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Abstract 

The main purpose of this paper is to study the structure of well-known non-iterative Multiple Signal 
Classification (MUSIC) algorithm. We investigate a relationship between MUSIC-type function for imaging 
of extended thin electromagnetic inclusions or perfectly conducting cracks and Bessel function of integer order 
of the first kind. This is based on the structure of left-singular vectors of collected Multi-Static Response 
(MSR) matrix whose elements are measured far-field pattern. Some numerical examples are illustrated for 
supporting our investigation. 
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(-H ■ electromagnetic inclusions, perfectly conducting cracks, numerical examples 

Cli 

1. Introduction 

One of the goal of inverse scattering problems is to identify unknown properties (e.g., shapes, material 
properties, locations, and constitutions) of targets from measured scattered field data. Generally, this 
problems are solved by the Newton-type iteration scheme; minimization of scattered data and computed 
data by generating an admissible cost functional. Related works can be found in P,.16, 22, 31] and references 
therein. Unfortunately, 1. non-convergence or appearance of several minima are arise in the reconstruction 
procedure due to the non-convex property of cost functional, 2. complex Frechet derivative must be evaluated 
at each iteration step, and 3. a priori information of targets is required. Although, above conditions are 
full-filled, reconstruction procedure will be fail if one begins iteration process with a bad initial guess. 
Hence, generating a good initial guess close to the unknown target must be performed preferentially. For 
this, alternative non-iterative imaging algorithms has been developed in order to obtain a good initial guss 



for ex amp le linear sampling method IJ, [21|, multi- frequency based Kirchhoff and subspace migrations 
0, H 111 H m, [ill , and topological derivative strategy [3, [lO, [Illli, [13 ■ 

In the pioneering research 15|, based on the eigenvalue structure of MSR matrix, MUSIC algorithm 
have been developed to find locations of point-like scatterers. Recently, MUSIC algorithm was applied to 
the various problems for example, detection of anti-personnel mines [5[, location search of small inclusions 
[6L[l3| and hidden corrosion in pip e line 0], shape reconstruction of arbitrary shaped thin inclusions, cracks, 
and extended targets 0,0,111) H, 2^- Based on these results, MUSIC algorithm offers an accurate location 
of small inclusions. In contrast, it is also effective in terms of shape identification but does not yield the 
complete shape of extended targets. Hence, obtained results adopted good initial guesses and successfully 
performed iterative based algorithm. 

Although, MUSIC algorithm offers good results for small and extended targets, a detailed structure 
analysis must be considered since some phenomena cannot be explained for example, unexpected appearance 
of ghost replicas along the boundary of targets [l^, Figure 9(b)] or an image with poor resolution [30l 
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Section 4.4]. The above phenomena gives us an motivation of exploring some properties of MUSIC-type 
based imaging algorithm in the shape identification of arbitrary shaped thin electromagnetic inclusions and 



cracks considered in [281 . 130j . Our analysis is based on the rigorously derives asymptotic expansion formula 



in the existence of thin inclusion and physical factorization of MSR matrix [19[. With this, we will 
establish a relationship between MUSIC-type imaging function and Bessel function of integer order of the 
first kind. 

This paper organized as follows. In section [21 we introduce two-dimensional direct scattering problems, 
asymptotic expansion formula in the presence of thin inclusions and MUSIC algorithm. In section [31 we 
identify structure of MUSIC imaging function and discuss some properties. Some numerical experiments 
are shown in section [4l This paper ends with a short conclusion in section [5l 

2. Direct scattering problems and MUSIC algorithm 

2.1. Two-dimensional direct scattering problems & asymptotic expansion formulae 

Suppose that an extended electromagnetic F with small thickness is located in a two-dimensional homo- 
geneous space We assume that F is characterized by its thickness h and supporting curve a as 

cr = {x -f 77n(x) : X e (7, t] e {~h, h)}. 

Throughout this paper, we assume that F and are classified by their dielectric permittivity and 
magnetic permeability at a given frequency uj. Let < Eq < -l-oo and < /xg < +oo denote the permittivity 
and permeability of K^, respectively, and analogously < e < -|-oo and < /i < +oo, those of the inclusion 
F. Then, we can define the piecewise constant dielectric permittivity e(x) and magnetic permeability /i(x) 
as _ _ 

s(x) = ( f - ^ f \r ^ r ,0 for X e E^F 

^ ^ [ e for X e F '^^ ^ [ /i for x e F, ^ ^ 

respectively. For convenience, we set Eq ^ — 1, e > Eq, and > ^q. 

For given fixed frequency lu (we assume that wave number k — lo^EqUq = w), we denote by 

uo(x,0;cj) ~ e'"^-'', x e (2) 

a plane- wave incident field with the incident direction e where §^ denotes the unit circle. Let 
m(x, 0\ u)) — uo(x, 6, uj) + Us{'x., d; oj) denotes the time-harmonic total field which satisfies following Helmholtz 
equation 

V ■ Vm(x, 9; uj)^ + tj2e(x)M(x, 6>; w) = (3) 

and unknown scattered field Us(x, satisfies the Sommerfeld radiation condition 

lim Vrf — — ' — ^ ^ iu!Ush<.,0;uj) \ = uniformly in all directions ■; — r G 
r--!-o \ or ) ' |x| 

Notice that above radiation condition implies the asymptotic behavior 



Ws(x, 0; lS) = — ^=Uoo(x, 6;ijj) + o { — -= for all |x| — > +oo, x = 




X 



Then by virtue in 11| , the far- field pattern can be written as following asymptotic expansion formula 
«oo(y, e; uj) = /^^^^^ ^ ((-y) • M(x) • + (e - 1)^ e-(«-^)-da(x) + o{h). 



Here, a 2 x 2 symmetric matrix A(x) is defined as follows: for x e ct, let t(x) and n(x) denote unit tangent 
and normal vectors to a at x, respectively. Then 
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• M(x) has eigenvectors t(x) and n(x). 

• The eigenvalue corresponding to t(x) is 2 — 1 

• The eigenvalue corresponding to n(x) is 2 (1 — /i). 

2.2. MUSIC-type imaging algorithm 

Assume that we have the same incident and observation configuration. Then collected MSR matrix can 
be written 

where Kji is the far-field pattern in the j—th direction of observation and Z— th direction of incident: 



Kji = Uoo(yj,^/;' 



yr- 



4:\/ujn 



=h 



^2(1 + i) length(cr) 



4-v/w7r 



M 



(x)-0, + (£-l) e'"(^^+^')-"da(x) 



M 

E 



m— 1 ^ 

+ 2(1 - • n(x,„)6>, • n(x„) + (e - 1) 



2 ( i - 1 ) 6>j • t(x„06'i ■ t(x„0 



(icr(x), 



where length((T) denotes the length of a, refer to [30l |. 

The MUSIC algorithm is as follows. First, let us perform the Singular Value Decomposition (SVD) of ] 



M 



K = USV* « J2 5„,U„,(x™)V™(x„ 



where Um and V™ are left- and right-singular vectors of K and 5*^ denotes the non-zero singular values. 
Second, based on the structure of left-singular vector, we define a vector f(z) € C^^"'^ as 



f(z) 



ci • (1, 0l)^e-«-^ C2 • (1, 02)^e^"^-^ • • • , cjv • (1, 9Nfe 



iujO L -z 



(4) 



where the selection of c„ e {0}, n — 1,2, ■ ■ ■ ,N depends on the shape of the supporting curve (t(x). 
Finally, by defining a projection operator 



P._.(f(z)) := I 

we can design MUSIC-type filter functional: 

E(z) := 



UN 



M \ 

- ^ U,„U^, f(z), 

m=l ) 



|P._e(f(z) 



(5) 



Remark 2.1. Based on the results ]29i . \3ui]. the selection of Cn is a strong prerequisite. The selection depends 
on the shape of the supporting curve a, e.g., c„ = (1,0,0)^ is a good choice for a purely dielectric contrast 
case. But for a purely magnetic contrast case, c„ must be of the form c = (0,b)"^, where h is a linear 
combination of t(x„j) and n(x„i) for m = 1,2, ■ ■ ■ , M . Unfortunately, we have no a priori information of 
shape of a. Due to this reason, in f7^ . \2^]. a set of (large number of) directions are applied in order to find 
optimal vector b. Applying this b yields a good result but this process requires large computational costs. 
Hence, we assume that c satisfies c„ • (1,0„)^ ~ 1 for all n, i.e., 



f(z) = 



(6) 



and explore some properties of MUSIC-type imaging algorithm. 
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3. Properties of MUSIC-type imaging function 

In this section, we determine MUSIC-type imaging function. Before proceeding, we recall a useful result 
developed in [T^j- This plays a key roll of our explore. 

Lemma 3.1 {^). For sufficiently large N, e §\ and | € R"^ , 

1 ^ ^ r 

- ^ « - / e^-'--dS{e) = Jo(c.|x|), 

1 ^ ^ r 

5. J. T/im electromagnetic inclusion case 

For the first and simplest case, we consider the the dielectric permittivity contrast case, i.e., e ^ Eq and 
/i = /io- The result is follows. 

Theorem 3.2 (Only dielectric permittivity contrast case). For sufficiently large N and uj, can be written 
as follows: 

X -1/2 



1 



/ M 
V \ m— 1 



|P„„.(f(z)) 

where f — g means that there exist a constant 7 such that f = jg. 
Proof. In this case, the left singular vector is of the form 



(7) 



Notice that since oj is large enough 

N 



1 



_ 1 r 

U,„ • U,„, ^ ^ 5] g»-e„.(x„-x„,) ^ Jp(^|x^ _ x,„,|) = J 



1 if m = m' 
if m 7^ m'. 



(8) 



Then P^ais^ can be written as 

P„„..(f(z)) : 



HAT 



M \ 

- ^ U,„C f (z) 



-E 



M 



M 



^iiuO N --X-m ^iu6i-{z — 'X.^-n) _j_ g^LJ0 jv 'X^ g2a;02 ' ( z — X„i ) 



M 



^iujOjsi z 



m— 1 

since e*"^"'^ = e*"^"-^™e^"^" (^"'''"). Then 



ll^no..(f(z))|| = 



N 



M 



.n— 1 \ m— 1 



E 1- E e'"«"-^"""'"^+e-'"^"-^^-''"Mjo(^^|z-x„|) 



M 



J2 e^"^"-'''"Jo(a;|z-x„|) ^ e-^"^"-^- Jo(a;|z - x^|) 



1/2 



Since 



N M 



n—1 m—1 



M 



(9) 



m=l 



and applying orthogonal property ([5]) yields 



A/ 



A/ 



^ e^-«"-'" Jo(c^|z - x,„|) ^ e-^"«"-^'" Jo(w|z - x,„|) 

\m— 1 / \m — 1 

= ^ e*-«"-(-"'-™) Jo(c^|z - x^l) Jo(c.|z - x^-l) 

= A^Jo(w|x,„/ - x,„|)Jo(w|z - x„|)Jo(w|z - x„i'\) = 
Therefore, by (jH]) and (TlUl) . we can obtain following result 

||P„_(f(z))|H ViV (l - ^ Jo(..|z-x„|)M 

\ m=l / 



(10) 



A^Jo(w|z-x,„|)2 if 



if m ^ m'. 



1/2 



□ 



Note that Jo{x) has its maximum value 1 at x = 0. This means that ([7]) plots peaks of large magnitude 



(oo in theory) at the location x^ e a and small one at x ^ ct (see Figure 1(a) ). This is the reason that why 
MUSIC algorithm offers good result for purely dielectric contrast case for the full-view inverse scattering 
problem. We refer to Figure 2(a) and various results in 19, 28, 30]. 

Next, we consider the magnetic permeability contrast case, i.e., s — Sq and fi ^ /iq- The result is follows. 

Theorem 3.3 (Only magnetic permeability contrast case). For sufficiently large N and lo, 0) can be 
written as follows: 



1 



1 ( 



-1/2 



|i^„„..„(f(z))|| 

Proof. In this case, the left singular vector is of the form (see [3; 



(x„)))^ Ji(w|z - x,„|)^ 



(11) 



U 



2(m-l)+s 



01 • ^,(x„Oe'"''-''-,02 • |,(x™)e*-«--",. • • • C(x™)e 



Au9 jv -Xji 



where 



^^(Xm) 



x(x,„) if s = 1 
n(x,„) if s = 2. 



Then based on the orthonormal property of singular vectors, we can observe following: if uj is large enough, 
and if TO to' or s 7^ s" then 



U2(m-1)-Hs • ^2(m'~l)+s" — X/ ( ^" ' (^s(x™) + £,s" {^m') j 

n=l ^ ^ 



(12) 



i(x„j - X„') • (Cs(Xm) + ^5//(x„i'))Ji(w|x„i - X„'|) = 0. 
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At this moment, we select f(z) in ^ then P„oiac can be written as 
/ M 2 \ 



P.„.e(f(z)) 



^iujOi -z 



n=l 
N 



Af 2 



TV 



{On ■ |,(x™))e-««-'" 5](0„ • ^,(x„))e-«- 



(z-x„) 



M 2 



^ ^(01 • ^,(x„0)((z - x,„) • ^,(x„))e^-«--'" Ji(c.|z - x„ 



m— 1 s—1 
M 2 



e 



iuj92 -z 



- * 51 • C.(x„0)((z - x,„) • |,(x„0)e*"^^-"'" Ji(c^|z - x„,|) 



m— 1 s—1 



M 2 



E E(^^ • ^si^"^mz - x„0 • |,(x„0)e*-«"-'^" Ji(l.|z - x^l) 



m— 1 s — 1 



Then 



||P._e(f(z) 



AT 



J\/ 2 



E 1 + * E E(^" ■ ^«(^™))((^ - • l.(x™))e^"«-(^-'''") Ji(^|z - x™|) 



n—1 ^ m— 1 s—1 

M 2 



E E(^« • ^.(^™))((^ - '^"O • C.(x™))e-*'^^"-(^-^") Ji(w|z - x„,|) 



m— 1 s — 1 
Af 2 



E E(^« ■ - ^rn) ■ |,(x™))e-«"-- Ji(l.|z - x™|) 



m— 1 s — 1 
7\/ 2 



E E(^« • - x,^) . U^,r^))e-'-''""'-JiHz - x,„|) 



■ m— 1 s — 1 

Applying Theorem 13.11 we can obtain 

N , 



1/2 



E *(e„-^,(x,„))e-«"-(^-^-)-z(0„-^,(x„O)e-'"^"-(^-"-) = -2iV((z-x,„).|,(x„0) Ji(c.|z-x 



n=l 
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And, based on the orthogonal property (jl2[) . we assume that m = m! then 



({On ■ €.(x™))((z - x„0 • C,(x„0)e*-«" --"71(^.12 - x,,,!)") 



TV 



^-'"^"■"'"'Ji(^|z-x„,,|) 
■Y. ( • - x™) • |,(x„0)' + (0n • is"Mf{{z - x,„) • |,„(x„))2 ) Ji(a;|z -x„|)2 



n=l 



E(^" • U^rn)r^ - x„) . ^,(x))Vi(a;|z - x„|)2 



Ji(cj|z — x„ 



((z-x„0-^.(x))^ / (VP-|,(x„))2d5(<p) + ((z-x„0.|,„(x))2 /■ (vp.|,„(x„))2d5((p) 

(14) 



Now, we consider the polar coordinate; since e let <p = (cos 0, sin 0)'^ and — (cos -0, sin then 
elementary calculus yields 



{(p ■ ^,{yi,n)) dS{ip) = cos\(t) - ip)d(t) = 2 cos^ = 2 



Therefore, by p^ . p^ . and p^ . we can obtain 



isin0cos(/)+ i( 



(15) 



|P„_(f(z))|| = Va 1^1 - ^ ((z - x,„) • (t(x,„) + n(x„0))Vi(c^|z - x„|)2 j 



1/2 



□ 



Opposite to the case of permittivity contrast case, (jlip will plots two curves with large magnitude in 
the neighborhood of supporting curve a because Ji{x) has maximum value in some neighborhood of x = 0. 
This is the reason why two ghost replicas with large magnitude and many afterimages with small magnitude 
are appear in the map of (jlip instead of true shape of supporting curve a when we select c„ such that 
c„ • (1,0„)^ = 1 (see Figure 1(b)). Related numerical examples can be found in Figure 2(b) and in [l^ 
Section 5]. 

Finally, by combining ([7|) and ([TT]) . we can obtain following result 

Theorem 3.4 (Both permittivity and permeability contrast case). Let e ^ Eq and fi ^ fiQ. Then, for 
sufficiently large N and uj, (0) can be written as follows: 



ie,p(z) := ^1 - Y (^-^o(w|z - x„i|)^ + ((z - x,„) • (t(x„) + n(x,„)))^Ji(w|z - x,„|)^^^ 
This result shows that (|16p plots large magnitude at z satisfies z ^ x^, and 



-1/2 



(16) 



M 



Y JoH^ - x„|)2 + ((z - x„0 • (t(x™) + n(x™)))Vi(a;|z - f \ =1 



It means that a result with poor resolution will appear. See Figure 2(c) and examples in [30|, Section 4.4]. 

7 



3.2. Perfectly conducting crack case 

In this part, we let F be a smooth curve which describes the crack: for an injective piecewise smooth 
function : [-1,1] — ^ R^, 

r = {0(z):ze [-1,1]}. (17) 
Let u{x,6;uj) be the single-component electric field that satisfies the Helmholtz equation: 

Au(x,0;w) + w2u(x,0;a;) in R^XT. 

For the sound-soft arc (Transverse Magnetic polarization - TM) case, u(x, 0: uj) satisfies the Dirichlet bound- 
ary condition on F: 

u(x, 6;uj) — on F 

and for the sound-hard arc (Transverse Electric polarization - TE) case, u(x, satisfies the Neumann 
boundary condition on F: 

du{x,e;uj) 

— — -- — =0 on F, 
c'n(x) 

where n(x) is a unit normal vector to F at x. 

Then the far-field pattern Uoo(y,0;w) for the scattering of an incident field uo{x,6;uj) = e*"^ '' from F 
is given by (see [22i] ) 



l + i 

(i-»)V^ 



sound-soft arc 



y • n(x)e 6)dx : sound-hard arc. 



Based on the physical factorization in [l^, [2^ , for coincide incident and observation configuration, left- 
singular vector of MSR matrix is of the form: 



U„,, = 



(01 • n(x„0)e^"«i (02 • n(x„0)e*"«^-"", • • • , (0jv • n(x,„))e''^«"-"" 



sound-soft arc 



sound-hard arc. 



This means that the structure of left singular vectors for sound-soft and sound-hard arc is same as the one 
of permittivity contrast case and permeability contrast case, respectively. Hence, MUSIC imaging function 
is same as in Theorem 13.21 and Theorem 13.31 for sound-soft and sound-hard arc, respectively. 







1.25 






1.2 






X-axis 






1.05 







-1 1 

X-axis 




(a) Graph of (1 - Jo{uj\x\))~'^/'^ (b) Graph of (1 - Ji{u)\x\))-^/'^ 

Figure 1: Graphs of (1 - .Ji,{uj\x\))-^/'^ , v = 0,1, for uj = 27r/0.5. 
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Remark 3.5. Based on recent work i23l ]. the structure of so-called "subspace migration" is 

M 

Jo{uj\z — x„i|)^, permittivity contrast case 



IE.m(z) = < 



m—l 
M 



( (z — x-m) • (t(x„,) + n(x.,„)) Ji(w|z — Xm|) I , permeability contrast case. 

^ m— 1 ^ ' 

Hence, we can observe following relationship between MUSIC and subspace migration as 

1 



E(z) 



v/1 - Es«(z) 



if we select vector f (z) of (0) as 



f(z) 



T/iis result is exactly same as derived in Formula (7.4)]- 

3.3. Small electromagnetic inclusions case 

We briefly consider the MUSIC-type imaging for the imaging of electromagnetic inclusions S„i with small 
diameter r: 

Em x„j + r_B,„, to=1,2,---,M, 

where Bm is a simply connected smooth domain containing the origin. We assume that S„i are sufficiently 
separated from each other, and denote S be the collection of such inclusions. 

Same as section[21 we let u(x, B; w) satisfies ^ in the presence of S and uo(x, 0; w) is given by Then 
by virtue in [t'I, the far-field pattern can be written as following asymptotic expansion 



Uoo(y,0;a;) 



J2 \Bm\ i-y) ■ A(x™) -e + ie-l) e*-(«-y)-." + oir') 



(18) 



Here, A(xm) denotes the polarization tensor corresponding to the (see Q)- 
Then by a similar process, we can obtain following results: 

Theorem 3.6. Let N and uj are sufficiently large. Then (0j can be written as follows. 
1. Dielectric permittivity contrast case: 

-1/2 



1 / 

E,(z)~-= 1- ^ Jo(^|z-x^|f 

^ \ rn—1 



2. Magnetic permeability contrast case. 



M r n 2^ 

((x - Zm) ■ (ei + e2)) Ji(a;|z - x„|) 



m—l 



'1/2 



where {ei,e2} denotes an orthonormal basis ir, 



3. Both permittivity and permeability contrast case: 

M 



1 / 



Jo(cj|z - X,, 



((x - z,„) • (ei + 62))^ Ji(w|z - x,i 



-1/2 



Remark 3.7. Now, let denotes the perfectly conducting inclusions with small diameter. Then, based 
on the result in f71 . Theorem 3.1], asymptotic expansion formulas for TM- and TE-cases are similar to the 
Hence, structure of (0) will be same as Theorem[ 
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4. Numerical examples 



In this section, we illustrate some numerical examples. For this, we choose a thin inclusion F = {x + 
r]n(x) : X e (7, 77 G {—h, h)} with smooth supporting curve 



cr = |(z + 0.2,z3 + z2_o.3)^ : ze [-0.5,0.5]| 



For the illumination and observation directions 0„, they are chosen as 

/ 2mr 2mT^ 
a„= (^cos— ,sm — 

for n = 1, 2, • • • ,N. Total number of directions is A'^ = 40 and applied frequency is a; = 27r/A with wavelength 
A = 0.4. 

Figure [2] shows the map of E£(z), Ep(z), and Ee^^(z) in the existence of F. Based on this result, we 
can easily observe that MUSIC algorithm offers very accurate result for permittivity contrast case. For 
permeability contrast case, as we explored in Theorem 13.31 we cannot obtain true shape of F. However, we 
can obtain an approximate shape of F based on the identified two ghost replicas since F lies between them 
based on Theorem 13. 31 For both permittivity and permeability contrast case, we can recognize the existence 
of F but it is very hard to identify complete shape. 





3000 
2500 
2000 
1500 
1000 
500 



(a) Map of E£(z) 



(b) Map of E^(z) 



-1 -0.5 0.5 1 
X-axis 

(c) Map of E£,mu(z) 



Figure 2: Image result via MUSIC algorithm. White colored line is true inclusion. 



For numerical example of perfectly conducting crack, we select a smooth curve of (jl7|) as 
cp[z) ~ < [ z,- COS — H — sm — cos — — : z e [— 1, IJ > 




and same test configuration as above. Maps of E(z) are illustrated in Figure [3] By comparing results in 
Figure m we can observe Theorem 13.21 and 13.31 holds for sound-soft and sound-hard arc case, respectively. 



5. Concluding remarks 

Based on the structure of left-singular vectors of MSR matrix, we investigated the structure of MUSIC- 
type imaging function by linking a relationship with the Bessel function of integer order of the first kind. 

Recent work [25] has shown that the MUSIC-type imaging algorithm cannot be applied to limited-view 
inverse scattering problems. In order to analyze the structure, the integration in Lemma [3.1l on the subset of a 
unit circle must be evaluated; however, this evaluation is very difficult to perform and in our knowledge, there 
is no finite representation about it. Therefore, extension to the limited-view inverse scattering problem will 
be an interesting research topic. Moreover, the considered method can be extended to a three-dimensional 
problem and we believe that the MUSIC-type imaging function will be linked to the spherical Bessel function. 
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